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Abstract 

We have investigated skewed parton distributions in coordinate space. We found 
that their evolution can be described in a simple manner in terms of non-local, 
conformal operators introduced by Balitsky and Braun. The resulting formula is 
given by a Neumann series expansion. Its structure resembles, for all values of 
the asymmetry parameter, the well-known solution of the ERBL equation in the 
momentum space. Performing Fourier transformation we have reproduced known 
results for evolution of momentum-space distributions. 



*work supported in part by BMBF 
t Alexander von Humboldt fellow 



1 Introduction 



Recently, a lot of effort has been invested in exploring deeply- virtual Compton scattering 
and hard, exclusive meson production processes, see []TJ, and @, §|, §] for recent 
reviews of theoretical and experimental situation. The celebrated factorization theorems 
@, H, Hi assure that when photon virtuality is large the underlying photon-parton sub- 
processes are indeed dominated by short distances and can be calculated perturbatively. 
The necessary information about long-distance, non-perturbative dynamics enters in the 
form of distribution amplitudes of the produces mesons, and new, generalized ]I0, [8] 
parton distribution functions of the hadronic target. Factorizable hard, exclusive processes 
provide therefore a possibility to study new aspects of nucleon structure, which cannot 
be accessed in inclusive measurements. 

More detailed studies of hard, exclusive meson production [T2|, [HI], [T3|, |T^, |TB[ have 
demonstrated that cross-section of these processes should be sufficiently large, so they can 
be investigated in the next round of dedicated experiments. These estimates are based 
on leading order (LO) QCD calculations which typically neglect the scale-dependence 
of generalized parton distributions, except for a diffractive region Hl7 , 18| . There, one 
can argue [12, IS] that generalized parton distributions which determine the dominant, 



imaginary part of the amplitude are proportional to the corresponding forward parton 
distributions. On the other hand, kinematical conditions for future measurements are 
such that a lot of data will be collected outside the diffractive region. Once data points 
achieve sufficient accuracy, a practical discussion of scale dependence of amplitudes of 
hard, exclusive processes will become unavoidable. Evolution equations for generalized 
parton distributions can be found in pi |2T|, 



r0|] . Their general properties have been 
discussed at length in the literature |10|, KL p2l, and by now various numerical algorithms 



23| , |24| , 0, [25|, and theoretical methods |]27| have been proposed and tested numerically. 



Most of these methods make more or less explicit use of the local conformal operators 
which are multiplicatively renormalized at the one-loop level [28, 29]. Note also that so 



far the main effort has been devoted to study evolution of skewed parton distributions 
themselves, while away from the diffractive region the real issue is the scale dependence 
of the QCD amplitude. 

The goal of this paper is twofold. First, we aim at a new, practically applicable 
method of treating the problem of scale-dependence of skewed parton distributions and 
corresponding amplitudes of hard, exclusive processes. Second, we want to explore the 
explicit operator solution to one-loop QCD evolution equations found by Balitsky and 
Braun, which involves non-local conformal operators. Although it has been found about 
a decade ago, we are not aware about any previous attempt to make phenomenologi- 
cal use of it. As the Balitsky-Braun solution was formulated in terms of coordinate-space 
quark string operators it is natural to introduce here a concept of coordinate-space skewed 
parton distributions which are generalization of coordinate-space forward parton distri- 
butions [30]. As we shall argue below, hard exclusive processes involve momentum- and 
coordinate-space parton distributions in a rather symmetric way, so both languages are 
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equally convenient from the point of view of practical calculations. 

It turns out that the language of coordinate-space string operators can indeed be 
applied successfully to the problem of finding a solution of the evolution equations for 
coordinate-space distributions. The result can be interpreted as a Neumann series expan- 
sion of the latter in terms of matrix elements of multiplicatively renormalizable, non-local 
conformal operators. It can be understood also as an operator identity - an expansion of 
QCD string operator in terms of non-local conformal operators. 

Having found a solution of the evolution equations in coordinate-space, we have inves- 
tigated also a corresponding scheme for evolution of momentum-space distributions. Here, 
we have found that coordinate-space considerations suggest splitting of initial conditions 
into two pieces, corresponding to initial parton distributions with average momentum 
larger and smaller than the asymmetry parameter. Because of the linear character of 
the evolution equations, both terms evolve independently. The resulting algorithm is a 
combination of expansion in terms of local conformal operators of one part, and general 
expansion in terms of orthogonal polynomials of the other. Both techniques have been 
already discussed in the literature |TT], |T^|, Here we contribute to this discussion 
by showing how the corresponding formulae obey general self-consistency constraints. 

Our presentation is organized as follows. In the following section we introduce a 
basic framework for coordinate-space skewed parton distributions. In the next section we 
present a solution to the problem of scale dependence of the coordinate-space distributions 
as an expansion in terms of Bessel functions. Armed with this result, we show that it 
can be interpreted as an expansion in terms of matrix elements of non-local, conformal 
operators and deduce the corresponding operator identity. Our result is closely related 
to the explicit operator solution found in [p0[| . The equivalence of both approaches is 
explicitely demonstrated in the Appendix A. 

In the next section we turn our attention to momentum-space quark distributions. We 
will argue that coordinate-space considerations suggest a particular algorithm for treating 
evolution in the momentum space, and discuss general constraints which any such solution 
has to fulfill. The next section contains a summary. Finally, two appendices are devoted 
to mathematical details which are not discussed in the main body of the paper. 



2 Coordinate-Space Skewed Parton Distributions 



In this section, we introduce the basic theoretical description of coordinate-space skewed 
parton distributions. Consider the twist-2, gauge invariant, light-cone quark string oper- 
ator normalized at a scale /i 2 : 



0(a,f3) = q( — ~^-^-z)z 



a + (3 a — f3 
z, z 



z) z 2 =0 



(1) 



Here, [a, b] denotes the path-ordered exponential 

[a,b] = Vexp[—ig / A fl (x)dx^ 

Jb 
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which reduces to 1 in the Schwinger gauge z ■ A = (g stands for the strong coupling 
constant and denotes the gluon field). Coordinates a and (3 are defined in such a 
way that a z defines the center of the string composed from quark fields and the gluon 
line between them while (3 z corresponds to its 'length', understood here simply as the 
difference between coordinates of the quark operators. 

Charge conjugation odd and even quark string operators are obtained from 0(a, (3) by 
taking its components symmetric, respectively antisymmetric in f3. Denoting the former 
and the latter by + (a,f3), respectively Q~(a,j3) one has therefore 



± (a, (3) = - (0(a, (3) ± 0(a, -(3)) . 



(2) 



Factorizable hard exclusive processes are described by matrix elements of 0(a,/3) 
between nucleon states (n(P',S')\ and \n(P, S)) with the corresponding four-momenta 
P, P' and spins S, S'. Here, the momentum space skewed parton distributions enter 
naturally fIQ| , . In the following, we have chosen to use convention introduced by Ji [^0 
In this notation, the matrix element (|l]) can be represented as: 



{n{P\S')\0{a,(3) \n(P,S)) 
N(P',S')zN{P,S)e- 1 



duF(u,!;; Li 2 )e 



iuP(P-z) 



+ 



N{P', S'] 



iM 



■N(P, S)e 



duK{u^-^ 2 )e iu ^ p - z) . 



-i 



(3) 



Here N(P, S) and N(P', S') denote initial and final nucleon spinors, respectively. The 
average nucleon momentum is denoted by P = {P + P')/"2, and the momentum transfer 
is r = P — P' . The asymmetry parameter £ is defined by r ■ z = 2£P • z. 

From now on, we will explicitly consider the skewed quark distribution F(u, £; /i 2 ) only. 
Evolution equations for K(u, £;/i 2 ) are exactly the same, so one can replace F(w,£;/x 2 ) 
by K(u, £; ji 2 ) at any stage of the following presentation. To proceed further, we separate 
the symmetric and antisymmetric part of F(u, £; fi 2 ), such that 



iuP(P-z) 



duF(u,£,; fi 2 )e 



duF s {u^- fx 2 ) cos uf3{P-z) +i / du F A (u, f; /j 2 ) sin u(3{P ■ z) (4) 
L J jo L J 

Obviously, matrix elements of + {a,f3) are parametrized by F s (u, £; /i 2 ), or 'valence' 
(quark minus antiquark) skewed quark distributions, while matrix elements of 0~(a,/3) 
are given by F A (u, £; /i 2 ), or 'quark plus antiquark' distributions. 

Note that due to charge conjugation and time-reversal invariance, F s (u,l;; fi 2 ) and 
F A (u, £; fi 2 ) contribute to the corresponding amplitudes M 5 (£; Li 2 ) and M A (£; Li 2 ) of hard, 
exclusive processes through 



M 5 (£;/i 2 ) ex f'duF^u,^^ 
Jo 



+ 



u — It 



IE 
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M A (£;/i 2 ) oc fduF^u,^^ 
Jo 



£ — u — ie £ + u — ie 



(5) 



Now we introduce, in analogy to the forward case, coordinate- space skewed quark 
distributions T S {P 1 £; /x 2 ) and J rA (P,£,; fi 2 ) as 



^ 5 (/U;/i 2 ) 

The inverse transformation reads 



2 /-i 



7T JO 

2 



Jo 

[ dojF A (u,^n 2 ) sinM). 
Jo 



(6) 



dpF s (P,&n 2 ) cos (up), 

o 

dpF A {p,£-^ 2 ) sin (u/3) . 



(7) 



Inserting these formulae into (|5|) one finds that the amplitudes M s (£; /x 2 ) and M (£; /x 2 ) 
can be calculated directly in terms of coordinate-space skewed quark distributions as 

M 5 (£;/i 2 ) oc ct/ ^e-^^ 5 (/3,£;/i 2 ) 
Jo 

poo 

M A (£;/i 2 ) oc ztt/ dpe~ iK T A {p^[i 2 ). (8) 
Jo 

Note that everywhere above we have used the notation P for a variable which arises 
from the string length P, introduced in ([!]), by a rescaling /3 — > P(P ■ z). 

Equations @ and @ show that in general momentum- and coordinate-space distri- 
butions play a symmetric role the in description of hard, exclusive processes. This is a 
new situation; in the case of DIS only the former have a direct physical interpretation in 
terms of observables. In particular, equations (H) can be reinterpreted as 

roo 

M s {i-^ 2 ) oc (n(P',S')\ dPe-^ p - z O + (0,P)\n(P,S)) 

Jo 

POO 

M A (£;/, 2 ) oc (n(P',S')\ dPe-^ Pz O-(0,P)\n(P,S)) (9) 

Jo 

We have found it interesting to note that it is as convenient to discuss hard, exclusive 
processes in terms of the correlation function (Q) at various longitudinal distances, as it 
is in terms of momentum-space parton distribution functions. 

For numerical calculations it is often more convenient to rewrite the amplitude M A (£; /i 2 ) 

as 



M A (£; /i 2 ) oc \ Cdu u F A (u, f ; /i 2 
4 Jo 



+ 



£ — u — ie £ + u — ie 
The corresponding expression in terms of coordinate-space distribution reads: 



roo 

M(f; fx 2 ) A oc m \ dp e~ iK F A '(B& $ 
Jo 



(10) 



where we have introduced a special notation 

It is easy to see from (|j) and (0) that coordinate-space skewed distributions introduced 
here are generalizations of coordinate-space parton distributions known from the forward 
case. On the other hand, the physical interpretation of the coordinate-space skewed quark 
parton distributions is twofold. First, they describe, according to (||), non-forward matrix 
elements of the non-local string operator Q(a,{3). Second, through they directly 
determine amplitudes M s (^]fj, 2 ) and M A (£; (i 2 ). Comparing with the forward case one 
finds that the way nucleon structure can be probed in hard exclusive processes is indeed 
complementary to what happens in deep-inelastic scattering experiments, as the latter 
are sensitive only to the imaginary part of the amplitude M A Q We observe that real and 
imaginary parts of amplitudes (§) indeed probe in different ways, respectively through 
cosine and sine Fourier transformations, twist- 2 quark correlation function in a nucleon 
target at various longitudinal distances. 



3 Evolution of coordinate-space skewed parton dis- 
tributions 

In this section we discuss how the explicit scale dependence of coordinate-space skewed 
parton distributions can be derived using arguments based on the conformal symmetry 
29|j . Here we consider flavour non-singlet quark distributions only. A generalization to 



the flavour-singlet case, including effects of mixing between quark and gluon operators in 
the flavor- singlet sector, will be given in a separate paper. 

Explicit solutions of evolution equations for momentum-space skewed parton distri- 
butions rely on the fact that, due to the conformal symmetry, one can find basis of local 
operators which are multiplicatively renormalized at the one-loop level. Our aim here is 
to construct similar expansion in the coordinate space, but in terms of matrix elements 
of multiplicatively renormalizable non-local string operators. As we shall see, the corre- 
sponding coefficient functions are given by Bessel functions J M (x), and the whole procedure 
can be understood mathematically as the Neumann series expansion of coordinate-space 
distributions. 

Recall that a function f(x) can be expanded in a Neumann series according to [3^] 

00 r°° dX 

/(x) = £(2z/ + 2 + 4n)J, +1+2n (x) / — /(A)J„ +1+an (A). (12) 



1 In the forward, DIS case, the skewed quark distributions become the usual quark densities i.e., 
F A (u, M 2 ) -> F A (u, 0; y 2 ) = q(u; M 2 ) 
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In particular, one finds that {j3^) v e luJl3 can be decomposed according to the Sonine's for- 
mula |32||: 



{PZYe** = rT M E *> + ") C n^/0 Ju+niPO (13) 

n=0 

where T is the Euler gamma function and denotes the Gegenbauer polynomial, respec- 
tively. Choosing v = 3/2, and taking the real part of ( |TB"D one obtains: 



2 



= _ r 



oo 



J2 (-l)»(3/2 + 2n) ^ 2 J 3 / 2+ 2n(^) (14) 

n=0 



Inserting this expansion in the definition of the coordinate-space skewed quark distribution 
(0) and interchanging summation and integration one obtains: 



3 

^(/3,£ ;/i 2 ) = £(- 1 ) n (3/2 + 2n)J (/?£) / duF s (u,^^)C 3 2 i\u/0- 

VP?/ n=0 J ° 

(15) 

Now, note that a Gegenbauer moment || 

^ f X duF s {u^^)ClL 2 {u/^) 
Jo 

is proportional to the matrix element of multiplicatively renormalizable local conformal 
operator [^, |29| . Its scale dependence is therefore given by 

i 2n fdu F s (u,^; Q 2 ) Cln(u)/£) = fdu F s {u^^ 2 )C ? ^{u/i) (16) 

JO Jo 



where 

/ logtQVA 2 ) ^ 
" llog^VA 2 )^ 
Here, A is the QCD scale, ^{k) is the anomalous dimension, 



(17) 



7(*) = 4/3 + j^—^ - 4(¥(* + 1) + 7 e)) , 

with ^(sc) = dL °^ ; 7b being the Euler constant, and b = 11 — 2/3Np. 

As it follows, the scale dependence of the coordinate space distribution J 78 ^,^; fi 2 ) is 
given simply by 



1 / 2 \ 2 00 /"l 

^ 5 (Ae;Q 2 )=-^ (^j E^ 1 )"^^^^) y o ^F 5 (u;,e;/i 2 )c 2 3 / 2 (^). 

(18) 



6 



In the case of T A {f3 1 ^ n 2 ) one can apply a decomposition of sin (uf3) analogous to 
141): 



2 



3 



- I 



3 



oo 



J2 ("l) n (5/2 + 2n) C'iUuj/O J 5/2+2n m (19) 

n=0 



to obtain 



i / o \ 2 00 rl 

^(/?,e;/i 2 ) = ^= f-^J E(- 1 ) n (V2 + 2n)J 5/2+2n (/30 j / o ^F A (^,e;/i 2 )C 2 3 / 2 +1 (^/0. 

(20) 

If JF A (/?,^; /i 2 ) is associated with a flavor non-singlet quark distribution, 

Jo 

is again multiplicatively renormalizable, but instead of ([IB]) we have || 

e 2n+1 jfdw F a (lo, £ Q 2 ) Cli 2 +1 (u/0 = L 2n+2 ? n+1 f/cu F A (u, £; /, 2 ) C^w/O . (21) 
The resulting scale dependence of JF A (/3, ^; Q 2 ) is governed simply by 



^(/UiQ 2 ) -57 E(- 1 ) n (V2+2n)L 2n+2 J 5/2+2n (/?0 / dw £; » 2 ) C 2 3 f +1 (o;/£) 

(22) 

Now we show that equations ( [Dp and (|22|), can be naturally understood as expan- 
sions of coordinate-space skewed quark distributions in terms of matrix elements of non- 
local, multiplicatively renormalizable, conformal operators. Indeed, applying flT2"|) one can 
rewrite fll8|) and (p2|) Neumann-type series: 

F S (P,Z;Q 2 ) = /3-^£(3 + 4n)J 3/2+2ri (/3e)Wi I™ d\yf\F s (\i-ii 2 ) J 3/2+2n (\£) 

n=0 J ° 

OO n QQ 

F A {P,Z;Q 2 ) = P~ l2 J2( 5 + 4n ) J 5/2+2n(f30L 2 n + 2 d\^\ f A (X, £; fl 2 ) J 5/2+2n (XO 

n=0 J ° 



(23) 



We are now in position to make relation to conformal symmetry explicit. Let us start 
from the obvious identity 



0(a,P)= da' dp'5(a-a')6(P-p')0(a',p'). (24) 

J-oo JO 

Applying (JT3) one finds a representation of a 5-function in terms of a Neumann series 

oo 

P6(P -/?') = E( 2 ^ + 2 + 4n) Jv+i +2n (0)J„ + i + 2n(0') ■ (25) 



ra=0 



Inserting this expansion into (|24| ) one finds that the string operator 0(a,j3) can be de- 
composed as 



0(a,/?) = /T 5 / -e- ifcQ E(l + 2j)Ji/2 + ,(|fc|/3)^(l/2 + j,A;; / i 2 ) (26) 



dk 

n=0 

in terms of conformal string operators S(l/2+j, k; /i 2 ). Here we have introduced a variable 
j = 1/2 + v + 2n which will be identified below with the conformal spin. Operators 



00 

ika 



S(l/2 + j,k;^) = dae lka dj3 J f3J 1/2+j (\k\(3) 0(a, (3) , (27) 

J — oo Jo 

introduced first in 1 201 , form a representation of a conformal group and are therefore 
multiplicatively renormalizable at a one-loop level |[20|| , i.e. 

5(1/2 + j, k- Q 2 ) = Lj 5(1/2 + j, k- /i 2 ) . (28) 

Because of symmetry properties in the variable /3, one finds that for charge-conjugation 
odd and even components of 0(a,/3), see equation (H), v has to be equal to 1/2 and 3/2, 
respectively. Equivalently, one finds that corresponding expansions run over operators 
S(l/2 + j, k; fj, 2 ) with the conformal spin j which assumes values j = 2n + 1, respectively 
j = 2n + 2. Combining this observation with (^BJ) and ( |2"5| ) one finds 

/oo AU „ 00 

— e- ika (3-i ^(3 + 4n)L 2n+1 J 3/2+2n (\k\(3)S(3/2 + 2n, k; /x 2 ) , 

-°° 271 n=0 

0-(a,(3;Q 2 ) = / — e- ifca /?-^(5 + 4n)L 2 „ +2 J 5/2+2 „(|A;|/?)5(5/2 + 2n,A;; /U 2 ). 

J-00 27T n=Q 

(29) 

Taking matrix elements of both sides of the above equations one immediately reproduces 
equations ( |TED and ([2"2]). 

Finally, let us note that equation ([23|) can be rewritten in the form 

^ 5 (A £; Q 2 ) = / rfA^ 5 (A, £; /i 2 ) £ L 2n+i(3 + 4n)* n 03, £)f n (A, (30) 

where = /^" 3/2 ^2n+i(/?0 and #„(A,£) = A 1/2 J 2n+ i(A£)- It corresponds to a 

standard representation of solution of a non-stationary evolution equation 

— F s (f3, £, r) = ^ d/?'#(/3, /3'; £)^ S (/3', £, r) . (31) 

Here, 'time' r is given by ^ log ( io|(^2/A2) ) • Functions \l/ n (/3,£) play a role of eigenfunc- 
tions of the corresponding stationary equation 

H^ n = E n ^ n . (32) 
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Functions \P n (/3,£) are eigenf unctions of the conjugated equation 

%H = Ejf n . (33) 
Energy levels are given by anomalous dimensions 7(2n + 1). The explicit form of the 



hamiltonian H can be found in []20| . Functions \& n and \l/ m are orthogonal with the scalar 
product given by the integral 

- r°° - S 

<* n |* m >=/ d/3* n (l3,0* m (l3,£) = -^ r . (34) 
Jo 3 + 4n 

The sum 

oo 

R S (P,P';£) = £L 2n+1 (3 + 4n)^(/U)*n(/3',£) (35) 

n=0 

is the standard representation of the resolvent of the evolution equation ([31]). Evolution 
of the charge-conjugation even quark distribution J rA (f3 ) £ t ;Q 2 ) can be interpreted in an 
analogous way. 

Note that coordinate-space conformal expansions (|23|) and ( P9f) have a form which 



closely resembles a general solution to the ERBL evolution equation [28, 33]. In mo- 
mentum space, however, such an expansion can be written only for distribution ampli- 
tudes i.e., in the case £ = 1. In analogy with the latter, a phenomenological model for 
coordinate-space distributions and, ultimately, physical amplitudes can be obtained by 
choosing expansion coefficients 



POO 

a s n (C,fi 2 ) = / d\V\F s (\,£;fi 2 )J 3/2+2n (\0 

J 



3 

~ V " " r W 5 (c^;/i 2 )^V/£) 



2^Uj v ' Jo 

POO 

«n(£,/i 2 ) = J o d\V\F A {\ £j /i 2 ) </5/2+2n(A£) 

, „ 3 

1 /2\ 5 , _ /■! 



2v^ U 




dwF^Cw.ej^cSS-iCw/O (36) 



according to some phenomenological model of low-energy nucleon structure. In partic- 
ular, it would be interesting to try a model in which (X n S £1X6 constructed from linear 
combinations of their values corresponding to the forward, £ — > 0, and totally exclusive, 
£ — > 1, limits. The resulting momentum-space distributions can be obtained by applying 
formalism developed in the next section. 

Equations (p3[) and ( f29]) , which provide a natural solution to the problem of LO scale- 
dependence of coordinate-space skewed quark distributions, belong to the main results of 
the present paper. Their form corresponds to an expansion in terms of a set of orthogonal 
eigenfunctions which are solutions of the stationary equations (j32|) and (|33|) . Note that 
a different representation of the same solution has been found in ||20||. The latter should 
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be understood as a representation of the Neumann series by an integral over complex 
values of conformal spin j. In the Appendix A we shall show in details how the Neumann 
expansion fl23|) can be obtained starting from the solution found in |2(| . 

Although we intend to discuss numerical aspects of the above algorithm in a separate 
publication, a short comment on its practical applicability is in order here. We have 

250 h 



15C 



im 5 (oi 2 



50 




Figure 1. Typical results of evolution of |M S (£)| 2 as a function of £, starting from a £- 
independent initial conditions F s (u,£; fift) = 1.1641iz -1 / 2 (l — u) 3 - 5 . The solid line denotes 
\M (£)| 2 at the initial scale fio = 1.777 GeV, the dashed line represents \M S (£)\ 2 evolved 
to fi = 10 GeV. 

checked, using various models of skewed quark distributions, that numerical algorithm 
for evaluation of physical amplitudes, based on equations ( p3|) and (H), gives accurate 
and stable results, see Figure 1 for an example, except for a case where the variable £ 
becomes small. This is related to an observation that a non-zero £ provides a natural 
cut-off for large j3 behavior of coordinate-space distributions, which significantly improves 
convergence of the Fourier integral (W) , as compared to the forward case pi . 



4 Momentum-space skewed quark distributions 

So far, we have discussed skewed parton distributions in coordinate space. In this section 
we turn our attention to the momentum-space distributions. Let us concentrate, for 
example, on F (u, f ; Q 2 ). It is related to F s {(3, £; Q 2 ) by the cosine Fourier transformation 
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(0). For simplicity of notation, we denote in this section F s (u,!;; Q 2 ) by F(u,£), and the 
corresponding distribution normalized at the initial scale /i 2 by F (u, £). 

Thanks to the conformal symmetry, the Gegenbauer moments of momentum-space 
skewed parton distributions are renormalized multiplicatively, but the corresponding Gegen- 
bauer polynomials do not form a complete set of functions on the interval —1 < u < 1 

A solution to the problem of scale dependence of F(u,£) based on an expansion in 
terms of orthogonal polynomials has been proposed in |f24] , [12|| and further developed, 
up to the NLO, in [p5| . Recently, another approach has been proposed in f27| . In the 
following, we shall discuss an alternative scheme, which originates from coordinate-space 
considerations presented in the previous section. 

In order to perform the Fourier transformation let us first split J 7 ((3, £) into two pieces, 
^<(/3,OaudJ>(/3,0: 

1 / 2 \ 3//2 00 f£ 
•F<(/U) = h7 E W(-l) fc (3/2 + 2fc)J 3/2+2fc (/30 jf dwFofoflC^w/fl 



1 / 2 \ 



fe=0 



3/2 



1 / 2 \ 

^((3,$) = — — ^L 2fc+1 (-l) fc (3/2 + 2fc)J 3/2+2fc (/30 / duF (u,OC 3 2 l 2 (u/0 

(37) 

Obviously, 

-F(/U) = ^<(/U)+^>(/U)- (38) 

It can be shown that both series expansions in fl3"T|) are convergent for all £ and /3. We 
denote their Fourier transforms as F<(m, £) and F > (u, £), such that 

F(«,0 = i ? <(«,0 + ^>(«.0- (39) 

At this stage, splitting of J 7 into JF< and .F> may seem arbitrary, but we shall demonstrate 
below that the associated decomposition of F into F < and F > offers new insights into the 
evolution of skewed quark distributions in momentum space. 

The next step is to insert the decomposition (|38|) into the cosine Fourier integral (0). 
In the case of F K (u, £) calculation is simple. Integrating corresponding series in fl3"T|) term 
by term one obtains the known result |H| 



(40) 

i duF^.OCl'^uj/i). 

The situation is, however, more complicated for ^F>(u,^) as the series obtained by 
performing Fourier integral diverges term by term. Instead, one can rewrite cos[w/3] as a 
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series expansion in Bessel functions, which follows from the Sonine's formula ([13]): 

cos[n/3] = cos[/3] + y/2rf (1 - u 2 ) £ (-^ ^^^^l^^) ^ 2 !^)^^^) • (41) 

The detailed derivation of this equation is given in the Appendix. Now, inserting the 
above expansion into (0) and interchanging the integration with summation one obtains 
a convergent series: 

oo n „j 

F>{u, = (1 ~ ^ 2 ) E JjS+t) ^ 2 («) E Wig(e a ) / dwFo^OCjfWO. (42) 

„ n 7. n *>c 



n=0 fc=0 ^ 

Here we have introduced the notation 



2 (-l)(^)r[3/2 + n + fc] 2fc 
fcl? J (n-fc)!r[3/2 + 2A;] 5 21 



n, 3/2 + n + fc 
5/2 + 2k 



(43) 



with 2-F1 being the hypergeometric function. It is easy to see that -R;£(£ 2 ) are polynomials 
in the variable £ 2 of order n for any k < n. 

As it follows, F(u,£) can be split naturally into two pieces, each of them being rep- 
resented by a series expansion in Gegenbauer polynomials. The first part, F K (u, £), is 
defined only in the region u < £. It has been noted |H| |8|, 24j that the evolution of 
F K {u, £) is governed by ERBL-type evolution equation p8|, |33|. The form of F<(u,£) is 
identical with a distribution amplitude 0(A) with A = it/£. The second part F > (u,^) 
can contribute for all < u < 1. This part is therefore similar to to the solution based 
on the expansion of F(u, £) in terms of polynomials orthogonal on a segment — 1 < u < 1 

Formula fl41|) can also be used to rewrite F<(w, £) as an expansion in terms of Gegen- 
bauer polynomials C^iu). In this case, the sum of F^u, £) and F > (u, £) reproduces 
solution found in 

Let us now demonstrate that F(tt, £), given by the sum of F<(w, £) and F > (u, £), and 
represented by series expansions fl40|) and (|42]), respectively, have all properties of the 
correct solution. For the sake of clarity we have omitted in the following all details of 
underlying calculations. Reader interested in mathematical aspects of this discussion will 
find the corresponding material in the Appendix B. 

1. Expansion of initial conditions. 

In the limit Q 2 = /i 2 , L 2 k+i = 1, one obtains F(u,£) = Fq(u,£). In this limit, two 
pieces in (^) give the following contributions: 

F < (u,£)Q2 =f t = 6{u<Z)F {u,Z), 

F>(«,0o>=^ = 0(u > OF (u,£) . (44) 

As Fq(u,£) is continuous and typically non-zero for u — £, it follows immediately 
that the function i*<(£, £) is nonzero despite the factor f 1 — ^-J in front of the 
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series in fl40|)! This means that, in general, the corresponding series expansion in 
Gegenbauer polynomials represents a function with the singularity at the end-point 
u = £. It belongs therefore to the functional space L p with 1 < p < 2. Further details 
related to the convergence of series expansion in terms of Gegenbauer polynomials 
can be found in [R^]. 

Note that the presence of the theta-function in fl40|) implies that partons which 
initially belonged to the segment < u < £ stay there in the course of the evolution. 
On the other hand, partons described by -F>(m, £), which belonged initially to the 
segment £ < u < 1, diffuse into the segment < u < £. An example of evolution of 
F > {u,£ t ) starting from ^-independent initial conditions is shown in Figure 2. It has 
been obtained taking 130 terms in the expansion ([42|). 

Moreover, considering evolution of Fq(u, £) to infinitesimally larger scale /i 2 — > /i 2 + 
<5/i 2 and interpreting the result as a new initial condition for the next infinitesimal 
step, one finds that upon evolution partons cross the border point u = £ from right 
to left side only - once they enter the segment < u < £, they never come back. 
These properties are in agreement with discussion in || [TTJ, based on explicit form 
of evolution equations. 




0.2 0.4 0.6 D.i.. 0.2 0.4 O.S 0.8 



Fig. 2a Fig. 2b 

Figure 2. Evolution of -F>(m,£ = 0.1) starting from a ^-independent initial condi- 
tions at the scale fi = 1.777 GeV, solid line. Figures 2a and 2b show F > (u, £ = 0.1) 
evolved to scales 10, respectively 1000 GeV, dashed line. Area under the dashed 
curve in the region u < £ represents number of partons which diffused from the 
region u > £. As the number of partons which initially occupied the region u > £ is 
conserved, this area is always equal to the area between dashed and solid lines for 
u > C 

Finally, both F < (u,C,) and -F>(m, £) can be treated as two independent solutions 
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corresponding to boundary conditions Hence, they should separately obey all 
constraints imposed by symmetry etc. on a valid solution. 

2. Asymptotic limit. 

In the limit Q 2 — > oo one finds: 

2 1/ v 2 \ r 1 
F(«,0 - ^<()-- 1-- ^F (u;,£). (45) 
q 2 ^oo 3 4 \ 4 2 / Jo 

In particular, the contribution of -F>(u,£) reads 

2 1/ ?y 2 \ r 1 
F > («,0 ^ e («<0g7 f 1 -^ / *"*b(w,£), (46) 



so in this limit all partons from the segment £ < u < 1 have diffused into the 
segment < u < £. 

3. Relation to solutions of ERBL and DGLAP evolution equations. 

In the limit £ — > 1 function F(u,^ = 1) is identical with the solution of the ERBL 
evolution equation for distribution amplitude: 

F(u, l)=F<(tt, l)=(l - u 2 ) £ L 2k+1 £ duF^ l)Cli\u) . 

(47) 

In the limit £ function -F(ii, £ = 0) is a solution of the DGLAP-evolution 
equation, represented as an expansion in terms of Gegenbauer polynomials C2^ 2 (w)Q: 

F{u , 0) = F>(u , 0) =(! - „*) ± (Jf^D ^C) x 
^ (n-Jfe)!(2Jfe)!r(3/2) Jo V ; 

(48) 

In particular one finds, as expected, that 

f 1 duu 2k F(u,0) = L 2k+l f 1 duu 2k F (u,0), F(u, 0) -* / duo F Q (u, 0) . 

(49) 

4. Conformal constraints . 

As required by conformal symmetry, Gegenbauer moments of the F(u,£) are mul- 
tiplicatively renormalizable || [2^, |29| : 

f 1 duF{u,0C 3 2 l\u/0 = L 2k+1 f l duF Q {u,i)Cli 2 (u/0- (50) 
Jo Jo 



2 This formula have also been obtained in Ei 
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As discussed above, the same is true separately for F K (u,£) and F > (u, £): 

f^duF^u^Cl^u/i) = L 2p+1 £duF Q {u>,i)C% 2 {Lo/ti) 
J^duF^OCl^u/O = L 2p+l £duF Q {u^)C% 2 {u)/i). (51) 

As L\ = 1, in the case of valence quarks the integral from F(u,£) over the whole 
domain is conserved and equal to the number of valence quarks in the target: 

f 1 duF(u, = f 1 duFoiu, f ) = N q . (52) 
Jo Jo 

In addition, the integrals involving _F<(tt, £) and F > (u,£) are separately conserved 
as well. 

5. Polynomiality. 

Moments of F(u,£), 

M n (0 = f 1 duu 2n F{u,0 (53) 
Jo 

can be related to matrix elements of usual, twist- 2 local operators of dimension 
3 + 2n. Combining this observation with time-reversal invariance and hermicity one 
finds that they have to be polynomials in £ 2 of a degree at most n. Note that this 
is a property of F(u,£), but not of F K (u, £) and F > (u,£) separately It is preserved 
by evolution. 

5 Summary 

The main goal of this paper has been to prepare a theoretical framework for a quantitative 
discussion of scale dependence of amplitudes of hard, exclusive processes with a general 
hadronic target. Factorizable amplitudes can be obtained either as a Mellin convolution, 
or as a Fourier integral, of momentum-space, respectively coordinate-space skewed parton 
distributions with corresponding Wilson coefficients. The explicit solution to the problem 
of the scale-dependence of coordinate-space skewed parton distributions, obtained in sec- 
tion |3], has a structure of a Neumann-series expansion of coordinate-space distributions 
in terms of matrix elements of non-local, conformal operators. Alternatively, it can be 
understood as the expansion of a QCD string operator in terms of non-local operators 
with a definite conformal spin. The particular form of both expansions explicitly reflect 
the conformal symmetry of one-loop QCD evolution equations. 

As discussed in [pU| , the evolution equations for skewed parton distributions in the 



coordinate-space representation have a simple form, analogous to the ERBL evolution 
equation for distribution amplitudes. Accordingly, the solution found here represents a 
series expansion in term of a set of orthogonal eigenf unctions of the evolution kernel in 
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the coordinate-space. It is related to the solution in terms of an integral over complex 



values of conformal spin j, found in pO 



Performing Fourier transformation between coordinate- and momentum-space parton 
distributions naturally leads to a scheme in which the momentum-space initial conditions 
are split into two parts, involving partons with average momentum fractions smaller and 
larger than the asymmetry parameter. The former never leave their initial domain, while 
the latter 'diffuse' in the course of the evolution from one domain to the other, speading 
over the whole interval. We have used our explicit solution to investigate this process 
numerically and found that the spreading occurs rather slowly. 

Preliminary numerical studies demonstrate that the coordinate-space approach leads 
to stable numerical results for scale-dependence of physical amplitudes away from the 
diffractive region. 
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Appendix A. 

In what follows we will demonstrate how the Neumann expansion (|T8|) of the charge- 
conjugation odd skewed quark distribution follows from the operator solution to one-loop 
QCD evolution equations found in [20|. As in Section |], we will denote in the following 
J 7 ([3, £; yU 2 ) simply by jF(/3, £). The authors of |2(J argued that conformal string operators 

dae ika / df3^Z 1/2+] (\k\f3)0(a,(3;n 2 ), (54) 

-oo JO 

where Z u is a solution of Bessel differential equation, form a representation of the collinear 
conformal group 5*0(2, 1) and are therefore renormcovariant. Moreover, applying the 
Lebedev-Kantorovich transformation Jjfj, they have shown that a flavour non-singlet 



QCD quark string operator 0(a, j3; Q 2 ) can be decomposed in terms of 5(| + j, k; /i 2 ) as: 

5+ica 

0(a, /?; Q 2 ) = -/T 3/2 / -£e~ ik <* J dj(j + |)J i+ i (1*1/3)^5(1/2 + j, k; /i 2 ) . (55) 



The Lebedev-Kantorovich transformation requires that S(h+j, k] (i 2 ), the multiplicatively 
renormalized conformal string operator (153), is defined by taking Z u = the Bessel 



function of the second kind p2[ . The contour of integration over j runs along the imagi- 
nary axis. 5 is a real number chosen in such a way that the contour lies inside the stripe 
parallel to the imaginary axis, where the integrand is an analytic function of j. 
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Matrix elements of the conformal string operator 5(1/2 +j, k; /i 2 ) can be easily found 
with the help of @: 



(n(P', S') \S(- + j, k; /i 2 ) \n(P, S)) = N(P', S') z N(P, S) x 



<2zr <)(k - T -.z) l n " d(3^(3Hl +j (k(3)(T s ((3,^fi 2 ) +tT A ((3,^fi 2 )) . 



(56) 



Finally, taking matrix elements of the both sides of the decomposition (55) we arrive at 
the following representation for the skewed quark distribution ,F(/3,£): 



5+ioo 



^,(;Q 2 ) = V 3/2 / rfiO' + |)^ +i (e/3)M(l+i,e;^), 



<5— ioo 



with 



S(i/2 + J ,e;/i 2 



d(3^(5'H 2 J+l/2 {i(3') / dcuF (cu,Z)cos[cuf3'}. 



(57) 



(5f 



To proceed further, we first integrate by parts over u. The boundary terms vanish because 
lim n ^ mFo(m, = -Fo(l,0 = 0. The resulting integrand has a regular large-/?' behavior, 
and therefore it is possible to interchange integrations over u and f3' to obtain: 



1 



7T 



"I) 



'2 



1/2 



r[i + jy2,i/2-j/2]/ (j,o+ 



i r[l/2+j,l/2-j/2,l/2-j] 



r[i + j/2] 



ii(i,0+ 



(59) 



,(_i )J+ 3/2 r[i/2 + j,i/2-j] . 

(l/2 + j)r[l/2-j/2] 2U ' U " 

To simplify the above expression we have introduced a special notation for functions 
I i(j ■>£,), i = 1,2,3, which are analytic in the region Rej > 0: 



/o(j,£) 
W,0 



l+j/2, 1/2 -j/2 
3/2 



a; 



r[i/2-i] 



2-^1 



1/2 -j/2, -j/2 
1/2 -j 



(60) 



^0(^,62^1 



l/+j'/2, 1/2 + J-/2 
3/2 + j 



^ 2 



Here Fq(u,£) = ^F (co>,£). The integration contour over j can be wrapped around the 
real positive axis, yielding the sum of residues in poles at j — 1 + 2k and j = 1/2 + k 
with k = 0, 1, 2, .... The series of poles at j — 1 + 2k originates from the first two terms 
in the rhs of ([5^). The second series of poles at j = 1/2 + k arises from the second and 
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third terms in the rhs of (|59|). We have found that residues of poles at j = 1/2 + k 
cancel exactly each other, so only poles at j = 1 + 2k give a non-zero contribution 
to the final result. The remaining calculations are tedious but simple. Finally, after 
integrating by parts over u in the opposite direction and rewriting the hypergeometric 
functions in terms of the Gegenbauer polynomials, the Neumann series flT8| ) arise. From a 
purely mathematical point of view this result illustrates a relation between the Lebedev- 
Kantorovich transformation and the Neumann series expansion, as discussed in ||36|| . 



Appendix B. 



Here we shall give the proof of equation (41) and provide necessary details of calculations 
related to discussion in section £|. 

In order to obtain formula (pEip we first expand cos(w/3) in the Neumann series using 
Sonine's formula (|T3|) with v = —1/2 and £ = 1: 

i R R °° i 

cos( M /3) = --(|)»r[-j] j_iC9) + (|)r[-i] E(-i) fe (-i/2 + 2k)c-J{u) J_i +2 M = 

k—l 

oo 

-1/2 



(61) 
one can 



cos(/?) + y/2n/3 £(-l) fe+1 (3/2 + 2k)C- { l%(u)J 3/2+2k ((3) . 



k=0 



rewrite C 2 (fc+i)( w ) as: 



Now, using known properties of the Gegenbauer and Legendre polynomials [35 

(-1) 



C2(fc+1)( M ) 



2(fc + 1 



-P 2 k+i(u) 



2(k + l 
(1-u 2 ) 



[uP 2k+ i(u] - P 2k (u)} 



d 



ClLUu) 



2(k + l)(2k + l)du 

— 1— (1 - u 2 ) Cl[ 2 (u) 

2(k + l)(2k+ 1) v } 2k \ ) 



2(k + l)(2k+ l)du 



(62) 



Inserting this formula into (|61|) one immediately obtains equation (41). 

Next, let us discuss properties of matrices (^ 2 ), introduced in equation (f43|). In our 
approach they originate from the cosine Fourier transformation of F > ((3,£), performed 
with the help of the decomposition (fUD. They can interpreted as expansion coefficients 
of a Gegenbauer polynomial C 2 p 2 (u) in terms of Gegenbauer polynomials C 2 ^{u/^) for 
arbitrary £, which for definiteness will be considered here to be greater than 0: 

4(«) = £^(£ 2 )4(«/£), 

fc=0 

m 2 ) = ,:,, L ;:w.:.".;,r ^ A 2/-j 



(p - fc)!T[3/2 + 2k] 



k-p, 3/2 +p + k 
5/2 + 2k 



(63) 
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Note that similar matrices appeared also in |53] as an element of a solution to the evolution 
equations for momentum-space skewed parton distributions based on the expansion of 
the whole skewed quark distribution in terms of polynomials orthogonal on the segment 
< u < 1. A straightforward way to obtain (|63"D is to take into account that polynomials 
C"2p 2 ( u /0 form an orthogonal basis for even functions of u on the segment — £ < u < £, 
with the weight function (1 — w 2 /£ 2 ). On the other hand, as C^iu) is a polynomial 
in u of a degree 2p, equation (|63"D ensures that left- and right-hand sides are identical 
independently whether u is smaller or larger then £. 

Note that the hypergeometric function ^ c l£ 2 ] with a = k — p, b = 3/2 + p + k 

and c = 5/2 + 2k is itself a polynomial in £ 2 of order p as long as a = k — p < 0, i.e. 
for all < k < p. In addition, the expansion inverse to ( |63"D is given by the same matrix 
function, but evaluated at l/£, i.e. R p (l; 2 ) has the following remarkable property: 

t R uve)RKe)=s lp . (64) 

k=l 

Unfortunately, we could not find a representation of R p k {^ 2 ) in terms of standard orthog- 
onal polynomials. 

Finally, let us concentrate on mathematical aspects which have not been discussed in 
details in section f|. 

Asymptotic limit: here one finds that as Q 2 — > oo, all L 2 k+i —>■ for k > 1 and L\ — 1. 
For F < (u,l;) only the first term in the expansion survives, and one obtains: 

2 1/ u 2 \ ft 

F<( W ,£) Q2=oo = 0(u<£)-- / o dw*b(",0- ( 65 ) 

Analogously, for F>(u,£) one obtains: 

F>KO Q2 =oo = (1 -u 2 )£<L,F (u;,Z) x g^^+M_^cJ?(«)i2ffK 2 ) . (66) 

3 € 



mine coefficients a p of the following expansion: 



To prove that the above series indeed equals to 6{u < £)I > (1 — u 2 /£ 2 ) one should deter- 



(ti < fl? 7 (1 - u 2 /e) = (1 - n 2 ) £ a n C 2 3 / 2 ( M ) (67) 

2 S n=0 



and compare them with the coefficients of the series in (|66|) . To this end, we multiply 

-y3/2 
'2/j 



both sides by Cf/ 2 ^) and integrate over u to obtain: 



(3/2 + 2p) 2 1 

a. 



--- J q du{l-u 2 /e)Cl>\u). (68) 



p (p+l)(2p + l)3£ 
The last integral can be easily calculated in terms of matrix i?j!(£ 2 ) 



£ ■ ^ £ 3 

f (1 - u 2 ie)Cg{u)dult = E *E(£ 2 ) / <V£ (1 - u 2 /e)Cl(u/0 = R p (e) , (69) 

k=0 Jo 
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so one indeed obtains 

a p = + W -rnf). (70) 

' (p + l)(2p+l) 



Expansion of initial conditions: when Q 2 = fi 2 , one should reproduce the initial distribu- 
tion F (u,£) as the evolution operator is equal to unity. Consider first F<(«, £). Setting 
L 2 k+i = 1 and interchanging summation and integration one obtains from ([401): 



W = *(« < 0~£du;F (u;,0 (l - ^ f) + ^ fc) g 2 f (VO^f (o>/g) ■ 

(71) 

Considering expansion analogous to ( |B"7| ) it is easy to see that the series above represents 
a 5-function: 

Ku/i - W /0 = (l " u 2 /e) p (J+^t + ^ cSfWOCS^WO , (72) 

and therefore one obtains 

F< (it, £)q 2 =^ = du FoK, £)% - m/0 = 0(u < F (u, f) . (73) 

Consider now F > (u,£),as given by fl42|) . In the case L2/C+1 = 1 the internal sum can be 
evaluated with the help of (|5"3|): 

*>(«,£ W = (1 - u 2 ) p 2 ^ 2n) Cli 2 {u) jT 1 dujF (u, flC^M = 

= jT 1 du F (u, £)6(u -u) = 6{u > F (it, 

Combining (|73|) and ( ff4|) one immediately finds the desired answer: 

F( U ,0q 2=m2 =F («,0- (75) 

Conformal constraints : here we check explicitly that conformal moments of F(u, £) 
are multiplicatively renormalized: 

f l duF{u^)C^ (u/0 = L 2k+1 f 1 duF Q {u,i)Cf (u/0 (76) 
Jo Jo 

The case of F < (u, £ ) is easy. Integrating fl42|) term by term one gets 

f duF^OCl^u/O = L 2p+1 fdcoF (cu,OC 3 2 i\uj/0 (77) 



(74) 
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(7S 



In case of F > (u, £ ) one can apply ( j63| ) and (|64| ) to obtain: 

f duF^OCl^u/i) = f>f(l/£ 2 ) f duF^OcTiu) = 
Jo l=Q Jo 

P n i 

= E<(ve 2 )Ewi5(e 2 ) / ^oK£)<? 2 ? (W0 = 

n=0 fc=0 ^ 

= e wi I 1 dwF { W ^)cii\w/oi:K{i/e)R n k {e) = 

k=0 ? n=fc 

= L 2p +i dwF (w, O c lp{w/i) . 

Sum of ([77]) and (|78|) gives explicitly (|76|). 

Polynomiality: Assume that moments M^(£) are polynomials in £ 2 of a degree n. 

Then, expanding u 2n in terms of C^p («/£) and using (|76|) one finds 

(79) 



fc^o k\T[2n — k + 5/2] 
The last integral can be expressed in terms of Mq(£) as 



e n f^F^Od^/i) = 1 g 2^-^r[2n -2k - p + 3/2] 2w _ 2fc _ 2p 
Jo ^ r[|J p\(2n - 2k - 2p)\ 

(80) 

Collecting all terms, one immediately finds that polynomiality is preserved by evolution. 

Finally, let us give some arguments in support of the convergence of the Gegenbauer 
series ([40]) and fl42]). Strictly speaking, we have shown that expansions are convergent 
for the case Q 2 = « 2 , or L 2 k+i = 1 only. We expect, however that evolution to higher 
scales will make the convergence only better because of the known fact that for large k, 
L>2k+i ~ k~ c with c > 0: 



L>2k+l ~ 

k— >oo 



const -IT^* MM. 
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